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Abstract 

The evolution of observable quantities of finite quantum systems is analyzed when 
the latter are subject to nondestructive measurements. The type and number of mea- 
surements characterize the level of decoherence produced in the system. A finite num- 
ber of instantaneous measurements leads to only a partial decoherence. But infinite 
number of such measurements yields complete decoherence and equilibration. Contin- 
uous measurements result in partial decoherence in finite time, but produce complete 
decoherence and equilibration as time tends to infinity. Resulting equilibrium states 
are characterized by representative statistical ensembles that, generally, retain infor- 
mation on initial conditions. Any system, to be observable, necessarily requires the 
presence of measurements, whose large number leads to the system equilibration and 
decoherence. 



PACS: 03.65.Yz, 03.65.Ta 



Keywords: Decoherence; Equilibration; Quasi-open systems; Nondestructive measure- 
ments 



E-mail address: yukalov@theor.jinr.ru 



1 



1 Introduction 



Decoherence and equilibration in quantum systems are the problems that has been studied 
in numerous papers. An extensive list of related literature can be found in the recent re- 
view articles [1-13]. Nowadays, this problem has gained much interest with regard to finite 
quantum systems. Such finite systems are now intensively studied both theoretically and 
experimentally because of their role in a variety of applications ranging from quantum elec- 
tronics to quantum information storage, processing, and computing [14-17]. Equilibration 
and decoherence from a strongly nonequilibrium initial state have been studied, e.g., for such 
quantum systems as spin assemblies [7,18-21], trapped atoms [6,13,22-27], and quantum dots 
[28]. 

One, generally, distinguishes two kinds of system dynamics: one is when the finite quan- 
tum system is isolated and another when it is connected to some environment. It is known 
that if the system is coupled to a sufficiently large equilibrium environment, it equilibrates 
and decoheres due to its interaction with surrounding [2,4,29-32]. But if the system is iso- 
lated, its motion is quasi-periodic and, thus, the system cannot equilibrate in the strict 
sense, displaying instead the Poincare recurrences [33]. However, an isolated system can 
equilibrate on average, relaxing to a quasi-equilibrium state defined by an ergodic average 
[34] and staying close to it most of the time [5]. Dynamics of isolated quantum systems de- 
pends on their closeness to integrability, revealing pre-equilibration effects [18,19,35,36] and 
displaying essential dependence on the way of their preparation, especially on the presence 
of defects [37,38]. 

In the present paper, we consider the intermediate case, when a finite quantum system is 
quasi-isolated, being almost isolated, except the action of nondestructive measurements, so 
that they do not destroy the system properties. We analyze the system decoherence under 
such nondestructive measurements. The main difference of the present consideration from 
the previous works is in the following. 

(i) The studied quantum system interacts not with an equilibrium bath, but with a 
nonequilibrium measuring device. 

(ii) The interaction part of the Hamiltonian is time-dependent, while a bath is usually 
described by a time-independent Hamiltonian. This essentially complicates mathematics 
and results in rather different consequences. The evolution operator now is not the standard 
exponential form exp(-Ht). Solving the evolution equation for this operator requires now 
to invoke the Lappo-Danilevsky theory. 

(iii) In the case of measurements, decoherence is not necessarily complete, as it would be 
for the bath, but the level of decoherence depends on the type and number of measurements. 

Specifically, nondestructive measurements are analyzed, but not arbitrary external per- 
turbations of the system. In that sense, the consideration is limited by exactly this type of 
measurements. Under this restriction, the overall treatment can be done for quite general 
conditions: the nature of the system can be arbitrary; it is not required that its spectrum 
be nondegenerate; the system states can be either pure or mixed; it is not required that 
the initial states be of the product type; no time averaging is needed; interactions of the 
system with the measuring device can be of arbitrary strength, but not necessarily weak; 
the measuring device can also be of rather arbitrary nature, provided it is nondestructive. 

As is stressed above, the main point is the consideration of the time-dependent interac- 
tion of the system with an external device, which makes the principal difference from the 
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standardly treated case of an equilibrium bath. Time-dependent interactions are typical of 
measurement procedures, when a measuring device is switched on and off. This is why it 
is possible to interpret such nonequilibrium interactions as measurements. Throughout the 
paper, the word "measurement" is used as a conditional term describing time-dependent 
influence on the system of an external nonequilibrium source. In addition, two limiting 
cases of the source influence, instantaneous and continuous, which will be analyzed in the 
paper, have the forms that are commonly associated with measurements [39-41]. Therefore 
the interpretation of the time-dependent source influence as measurement seems to be justi- 
fied. This, however, is not compulsory and one can treat the word "measurement" just as a 
brief term for the considered time- dependent interactions possessing several of the properties 
typical of measurement process. 



2 Quantum system under nondestructive measurements 

Let the quantum system of interest be described by a Hamiltonian Ha acting on a Hilbert 
space "Ha- The system is subject to a measurement procedure. The measuring device is 
characterized by a Hamiltonian H B acting on a Hilbert space % B . The total Hamiltonian is 
the sum 

Hab = H A + H B + H int , (1) 

where H int is the term describing the interaction between the system and the measuring 
device. Hamiltonian (1) is defined on the Hilbert space 

U ab = U a §§Ub. (2) 

Strictly speaking, the Hamiltonian (1) should be written as 

Hab = H A 0i A + iA0H B + H mt , (3) 

with 1a and 1 B being the unit operators on the corresponding spaces. However, it is com- 
monly accepted to omit the unit operators, writing, for simplicity, the Hamiltonian in form 

(!)• ' 

The measurement is called nondestructive when it does not disturb the system, in the 
sense that the system Hamiltonian is conserved, 

[H A , Hab] = . (4) 

Conversely, the studied system does not destroy the measuring device, so that its Hamiltonian 
is also conserved, 

[H B , Hab] = . (5) 
As a consequence of this definition, one has 

[Ha, Hi n t\ = [H B , H int ] = , [H in t, Hab] = ■ (6) 

The measuring procedure, enjoying these properties, can also be termed minimally disturbing 
measurement [42-44] or nondemolition measurement [45-48]. 
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The system Hamiltonian defines a complete orthonormal basis {\n)} by the eigenproblem 

H A \n) = E n \n) . (7) 
Respectively, the device Hamiltonian, by the eigenproblem 

H B \k)=p k \k), (8) 

defines a complete orthonormal basis {|&)}. 

Because of properties (4) and (5), the total Hamiltonian (1) possesses the eigenfunctions 
\nk) = \n) \ k). A measurement procedure is a nonequilibrium process, which leads to a 
complication coming from the fact that the interaction Hamiltonian Hi nt = Hi nt (t), generally, 
depends on time. Therefore the eigenvalues of the total Hamiltonian also depend on time, 
being given by the eigenproblem 

H AB (t)\nk) = [E n + e nk (t)} \ nk) . (9) 

The temporal evolution of the total statistical operator is prescribed by the rule 

PA B {t) = UAB(t)pAB(0)U% B (t) , (10) 

where the evolution operator is a unitary operator satisfying the Schrodinger equation 

i j UabM = H AB U AB (t) . (11) 

The total Hamiltonian H AB = H AB (t) depends on time through the interaction term 
Hi nt (t). This does not allow us to represent the evolution operator in a simple exponential 
form as it is usually accepted in the case of an equilibrium bath. However, we may notice 
that 



= 5 

Consequently, the Lappo-Danilevsky condition [49] 



H AB (t) [ H AB {t') dt' n 
Jo 

S kp [ E n + e nk {t) ] [ [E n + e nk (t') ] dt' . (12) 
Jo 



H AB (t), fH AB {t')dt' 
Jo 



(13) 



holds true, at least in the weak sense. Owing to this condition (13), one can represent the 
solution of Eq. (11) as 

U AB (t) = exp j-i jT H AB (t') dt'} . (14) 

The possibility of representing the evolution operator in form (14) is due to the Lappo- 
Danilevsky condition (13) that becomes valid when the measurement procedure is nonde- 
structive in the sense of conditions (4) and (5). 

Thus the total composite system is characterized by the statistical ensemble {H AB , p AB (t) } , 
with the statistical operator satisfying the normalization condition 

Tr ABPAB{t) = 1 , (15) 
where the trace is over the total space (2). 
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3 Temporal evolution of observable quantities 

What one is interested in any physical problem is the behavior of the system observable 
quantities that are represented by self-adjoint operators A composing the algebra of local 
observables A = {A} defined on the system space Ha- The measurable observables are given 
by the averages 

(A(t)) = Ti AB p AB (t)A , (16) 

in which the trace is over Hab- But, since A is given on Ha, the latter average (16) reduces 
to 

(A(t)) = Tr A p A (t)A , (17) 

with the partial statistical operator 

p A (t) = Tr B p A B(t) , (18) 

in which the degrees of freedom of He are traced out. 

Passing to a matrix representation transforms average (17) to 

W)> = EpIWA™ , (19) 

ran 

with A mn =< m\A\n > and a density matrix 

pi n (t) = (m\p A (t)\n) . (20) 

The latter can be written as 

pLW = 5>-,W, (21) 

k 

where 

pZk{t) = (mk\p AB {t)\nk) . (22) 

In general, any basis could be used for such a matrix representation [50]. But, as far as 
the system observables are of interest, it is convenient to employ the basis {\nk)} composed 
of the eigenvectors \n) of the system Hamiltonian H A and eigenvectors \k) of H B . Then, 
because of Eqs. (9) and (14), one has 

U AB(f)\nk) = exp |— iE n t — i J £ n k(t') dt'^j \ nk) . 

This, for matrix (22), gives 

Pinkif) = pinM exp j-iuW - I ^ Erankit') dt'^ , (23) 

with the initial condition 

PZM = (mk\pA B (0)\nk) , (24) 

transition frequencies 

Urnn = E m — E n , (25) 



5 



and the notation 

£mnk(t) = S m k{t) — £nk{t) ■ (26) 

By definitions (25) and (26), one has 

= , e nnk {t) = , (27) 
because of which the diagonal element 

PLit)=Y.PnnM = Pnn (28) 
k 

does not depend on time. Because of normalization (15), we have 

= £*» = 1 - (29) 

n n 

Separating in sum (19) the terms with m = n and m ^ n yields 

W)) = J2 P™ A ™ + £ PL(t)Anm ■ (30) 

As an initial expression Pab(O) one usually takes a disentangled factor product of Pa(0) 
and Pb(0). Here we do not assume this simplification, but keep the general form of Pab(0). 
Whether the latter is entangled or not does not play a principal role for what follows. 



4 Measurement procedure of several measurements 

In general, a measurement procedure can include multiple acts of measurement. Let there 
be M such measurement acts, so that the action of the measuring device be represented as 
the sum 

M 

in which fj(t) is a real function and Xj is a self- adjoint operator on Hab- 

According to the definition of nondestructive measurements in Eqs. (4) and (5) and its 
consequence (6), we have the commutators 

[H A Xj] = [H B , Xj] = . (32) 

Therefore the eigenproblem for Xj reads as 

Xj\nk) = ijnk\nk) , (33) 

with a real eigenvalue As a result, the eigenproblem for the interaction operator (31) 
takes the form 

H int (t)\nk) = a nk (t)\nk) , (34) 
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with the eigenvalue 

M 

<Xnk(t) = y~]€jnkfj(t) ■ (35) 
3=1 

In view of Eqs.(8) and (9), we get 



which leads to 
Then Eq. (26) gives 



£nk(t) = Ot nk (t) + Pk , 

£mk{t) — £nk{t) = Oi m k{t) — a nk (t) . 
M 

^mnfe(^) ^ ^ %jmnkfj(t) j (36) 
3=1 



■Ejmnk — Cjmk ^jnk ■ (37) 



where 

Introducing the notation 

Rmn(t) = Pmnfc(O) exp <j -i ^ X jmnk (pj(t) \ , (38) 

with 



M 



3=1 



<Pj(t) = I fj(t') dt> , (39) 
Jo 

for matrix (21), we obtain 

f£, n (t) = R mn (t) exp(-iuj mn t) . (40) 

5 Evolution after last measurement 

Function (39) characterizes the integral impact of the j-th measurement during the period of 
time [0, t\. Suppose that after the last M-th measurement, occurring at time t M , the integral 
impact (39) becomes 

= (t>t M ). (41) 

As is shown below, property (41) is valid for different types of measurements, including 
discrete measurements, whose action is equivalent to instantaneous kicking [51,52], as well 
as in the opposite case of continuous permanent measurements, acting uniformly in time 
[41,47]. 

Under condition (41), Eq. (38) reduces to 

Rmn{t) = P™M ex P {-ix mnk M<p(t)} , (42) 

k 

where 

M 

■Emnk — j^j- 

3=1 



1 M 
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Let us introduce the density function 

9mn{x) = Pmnk(W X ~ X mnk) • (44) 
k 

This function incorporates the properties of the measuring device, which affect the measured 
quantity. Therefore it can be called the effect density [39]. Definition (44) makes it possible 
to rewrite Eq. (42) as the integral transformation 

Rmnif) = J g mn (x) exp{-ixM(p(t)} dx . (45) 

The effect density (44) is normalized as 

j g mn {x)dx = Y,P™M. (46) 

J k 

Invoking the definition 

Pmn = Pl(0) = 5>^ fe (0) (47) 
k 

reduces normalization (46) to the form 

g mn {x) dx = p mn . (48) 
Therefore the effect density can be represented as 

9mn{.x) PmnPmn(,%) i (49) 

where the distribution p mn (x) is normalized to one, 

J Pmn(x) dx = 1 . (50) 

Thus for Eq.(45), we come to the expression 

Rmnif) = PmnD mn {f) , (51) 

with the decoherence factor 

D mn (t) = J p mn {x) exp{-ixMip(t)} dx . (52) 
Respectively, matrix (40) reads as 

Pinit) = Pmn(t)D mn (t) , (53) 

where 

Pmnif) = p mn exp(-iu mn t) . (54) 
As a result, for the evolution of observable quantities (30), we have 

(A(t)) = ^ ] pnnA nn + ^ ] Pmn(^)^4-nm-Dmn(^) • (55) 

This formula describes the evolution of observables after the system has been subject to a 
measurement consisting of a series of nondestructive measurement acts. 
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6 Decoherence caused by nondestructive measurements 



The second term in Eq. (55) is due to interference effects typical of coherent quantum sys- 
tems. Measurement procedure destroys coherence. In order to derive an explicit expression 
for the decoherence factor (52) one has to model the distribution function p mn (x) and to de- 
fine the type of the measurement procedure characterized by function ip(t). The measuring 
device is a macroscopic system, because of which its spectrum can be treated as contin- 
uous, similarly to the density of states of macroscopic statistical systems [53]. Therefore, 
the summation over k in the above formulas should be understood as integration over this 
multi-index. We shall consider two typical distributions, the Gaussian and Lorentz ones, 
and two opposite measurement procedures, instantaneous and continuous. This choice is 
based on the following arguments. A measuring device, being a macroscopic object, con- 
tains a large number of elements acting randomly on the system. As is known from the 
central limit theorem, the action of a large number of random elements is well described by 
Gaussian distribution. Because of this, the effect density of measuring devices is commonly 
represented by Gaussians [39,54-56]. Lorentzian distribution arises when measurement is 
realized by means of optical beams [40] . Two usually considered types of measurements are 
instantaneous [39,40] and continuous [39,41,55] measurements. 
In the case of the Gaussian distribution 



Here the standard deviation a, in general, can depend on the indices m and n. But for the 
simplicity of notation, this dependence is not shown explicitly. 
While for the Lorentz distribution 




(56) 



we have the decoherence factor 




(57) 




(58) 



the decoherence factor is 



D^(t)=e X p{-aM V (t)}. 



(59) 



When each act of the measurement procedure is instantaneous, such that 



/;(*) 



6(t 




(60) 



then Eq. (39) is the unit-step function 




(61) 



Consequently, after the last measurement at time t M , 



<p(t) = 1 



(t > t M ) . 



(62) 
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Then for the decoherence factor (57), in the case of the instantaneous measurements, we get 



a 



2 



D^ nst (t) = exp ( - — M'j , (63) 

where the indices m and n are omitted. And for the decoherence factor (59), we find 

Dt t (t)=e W (-aM) . (64) 

Another situation happens in the opposite case of a single but continuous measurement, 
when 

fj(t) = l, M=l, (65) 

so that 

<p j (t) = <p(t)=t. (66) 

Then the decoherence factor (57) is 

c? -' (t)=exp H(£> 2 }' (67) 

with the decoherence time 

tdec = ~ ■ (68) 

a 



And the decoherence factor (59) becomes 



D L cont {t) = exp (- -L) , 

V l dec J 



(69) 



with the decoherence time as in Eq. (68). 

In this way, a finite number of instantaneous measurements or continuous measurements 
during a finite interval of time lead to a partial decoherence in expression (55) of observ- 
able quantities. If the number of instantaneous measurements infinitely increases or the 
time of accomplishing continuous measurements tends to infinity, then decoherence becomes 
complete: 

lim ( A(t) ) = lim ( A(t) ) = V Pnn A nn . (70) 

M— >oo t->oo — ' 



Since the limit in Eqs. (70) does not depend on time, this also implies equilibration. 

The matrix p mn , according to Eqs. (22) and (28), generally, depends on the initial state 
of the total composite object, system plus device, since 

Pmn = ^2 (mk\ Pab(0) \nk) . 
k 

The dependence on the measuring device disappears if at the initial time the studied quantum 
system and the measuring device are not entangled, so that 

/W0) = Pa(0)6§p b (0) . 
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Then because of the normalization condition 

< k | p B (0) \k) = l, 

k 

the matrix element 

p mn = ( m | p A (0) I n ) 

contains information only on the system initial state. 

In the limit of infinite number of instantaneous measurements or large time of continuous 
measurement, the observable quantities reduce to the averages involving only the diagonal 
elements p nn , which demonstrates the importance of such diagonal terms [57]. The values of 
p nn can be found by defining the corresponding representative ensemble uniquely describing 
the system [58]. The diagonalization occurs not because of imposing some additional aver- 
aging over the random phases of initial states [59], or because of time- averaging resulting 
in ergodic averages [34], but happens naturally as a consequence of repeated measurement 
actions. 

7 Discussion 

The evolution of a quantum system, subject to the action of nondestructive measurements, 
is considered. Two types of measurement procedures are analyzed, instantaneous measure- 
ments and continuous measurements. A finite number of instantaneous measurements or 
continuous measurements during a finite time interval lead to partial decoherence. But if 
instantaneous measurements are repeated infinite number of times or a continuous measure- 
ment lasts infinite time, then there happens complete decoherence and equilibration. 

The interaction of a quantum system with a measuring device is principally different from 
its interaction with a bath. The latter is usually represented by a large equilibrium system. 
Contrary to this, the measurement procedure is a nonequilibrium process. Therefore the 
interaction part of the total Hamiltonian is time-dependent. This does not allow one to 
write down the evolution operator in the standard exponential form. To solve the evolution 
equation, we have to invoke the Lappo-Danilevsky theory [49]. 

We have considered the evolution of quantum systems under the action of nondestructive 
measurements. Then the system, starting form an initial nonequilibrium state, evolves into a 
state that is partially decoherent. The level of decoherence depends on the type and number 
of nondestructive measurements. Real measurements are usually supposed to be nondestruc- 
tive in order not to destroy the studied system. Of course, a measurement, being an external 
perturbation, could also be destructive. Moreover, it is possible to study the situation oppo- 
site to that treated in the present paper. It is admissible to investigate the problem of how 
an equilibrium system develops being subject to the action of a time-dependent external 
perturbation driving the system out of equilibrium [60]. But such destructive perturba- 
tions, generally, are not characterized as measurements. And this opposite process of system 
destabilization is a principally different problem. 

It is worth stressing that absolutely isolated systems do not exist. The notion of an 
absolute isolation is self-contradictory, since the statement that a system is isolated implicitly 
assumes that the fact of the system isolation is proved by observations. In the other case, such 
a statement cannot be true. To be true, a statement must be confirmed by observations. But 
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the latter, even being accomplished by nondestructive measurements, influence the system, 
leading to its, at least partial, decoherence. The absence of absolutely isolated systems and 
their decoherence, caused by observations, result in the system equilibration. 

The fact that the system isolation must be proved by measurements does not contradict 
the possibility of treating a system as being isolated during a finite time duration. For 
instance, one can prepare a system in a well defined initial state by a prescribed preparation 
method. Then, after a finite time, one accomplishes a measurement by a nondestructive 
state tomography and determines the final system state. Supposing that in between these 
two measurements, the initial and final one, the system has not been influenced by any 
other disturbances, it could be possible to assume that it has been isolated during this finite 
time interval. Such an assumption is nothing but a particular case of the finite number of 
measurements, as has been considered above. But, in any case, one has to accomplish at 
least two measurement procedures in order to assume that the system has been isolated for 
the finite time between these measurements. In such a case of a finite temporal observation, 
the system, as is explained above, does not equilibrate and exhibits only partial decoherence. 
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